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ABSTRACT 

A large fraction of bulgeless disk galaxies contain young compact stellar systems at their centers, in spite of the 
local gravitational stability of these disks. We evaluate two contrasting hypotheses for the origin of the nuclear 
star clusters in late-type disk galaxies. The clusters could not have migrated from distant eccentric locations 
in the disk. Instead they must have formed in situ, requiring radial transport of gas toward the center of the 
disk. This transport could be a consequence of the development of the magnetorotational instability in the 
differentially rotating warm neutral medium. We evaluate the rate of gas transport into the disk center and find 
that it is sufficient to support continuous star formation in that location. Enhanced stellar surface brightness 
in the inner few hundred parsecs and the formation of a compact stellar system in the central few parsecs are 
unavoidable in dark matter halos with divergent density profiles. We illustrate our conclusions on a model of 
the nearest late-type disk galaxy M33. 

Subject headings: galaxies: nuclei — galaxies: spiral — galaxies: star clusters 
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1. INTRODUCTION 



been growing in the compact 
in the centers of late-type spiral 



Recently, interest has 
stellar cl usters found 

galax ie s (ICarollo et al.l ll997t IBoker . van de r Marel & Vacc 
1999; Matthe ws et all Il999t IBoker etalJ 120011 1200:' 
Matthews & Gallagher 



. _ _, Boker, Stanek & van der Marel 

20031 ISchinnerer. Boker & Meierl 120031 IBoker et al J 120041) . 
Imaging with the Hubb le Space Tele scope of the central 
regions of ~ 75% ( Bo ker et al.l 12002) of late-type spirals 
(Hubble type Scd or later) has revealed young, compact 
stellar systems with luminosities -14 < Mj < -9 that dom- 
inate the photocenters of these disk galaxies possessing 
no discernable stellar bulge. The "nuclear star clusters" 
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(NSCs) have half-light radii R r } < 5 
masses - 10 4 - 10 7 M Q (Boker et al. 2002 
The dynamical masses of the NSCs, ~ 10 6 -10 !i M Q , 
an order of magnitude higher than would be expected from 
stellar synthesis models, indicating the presence of old er 
stellar generations that have faded (Walcher et al. 2004) 1 . 
Some of the NSC host galaxies exhibit spiral structure, dust 
lanes, and ot her morphological complexities, but many do 
not, leading IBoker et all (120021) to suggest that the NSC 
phenomenon is not related to nonax isymmetric features 
in the galactic disk. IBoker etafl ( 120041) find no evidence 
for a correlation between the presence of a cluster and the 
presence of a large-scale stellar bar. Stellar disks of the 
host galaxies have exponential surface brightness profiles to 
within hundreds of parsecs of the center. Surface brightness 
profiles of the inner few hundred parsecs but outside the NSC 
freq uently exhibit brightness excess abov e the exponential 
law (IBoker. Stanek & van der Marel2003l) . 

The formation mechanism of NSCs is unknown. The cen- 
tral rot ation curves of some Scd gal axies rise linearly with 
radius ( Matthews & Gallagher 2002), implying near solid- 
body rotation and a shallow potential well. Rotation curves 
of NSC host galaxies resemble those of low surface bright- 
ness and dwarf spiral galaxies. Radial mass profiles of the 
latter tw o classes of galaxi es have been exte nsively stud- 
ied (e.g.Jde Blok et al.l2003lls^ters et al.l2003USimon et al J 



2003). Models of the combined dark matter and baryonic d en- 
sity profiles that are finite at the center (e.g., Burk eril995l) . as 
well as those that are mildly divergent p oc r" 7 with 7 < 1 
(e.g., [Navarro, Frenk & White 119971) . both yield admissible 
fits to the rotation curves, while steeper profiles are ruled out. 

Two scenarios for the origin of NSCs are possible. They 
could have assembled in situ from the material available in 
the central region of the disk. Alternatively, they could have 
formed remotely and have subsequently migrated to the cen- 
ter via dynamical friction. We test these scenarios on a 
model (§ of the closest and the well-studied NSC host 
galaxy M33. In §|3] we argue that the first scenario — in situ 
formation — is the likely origin of NSCs. In §f4] we demon- 
strate that the alternative migratory scenario fails to deliver 
NSCs to their observed locations, and we discuss observable 
implications of our model. 

2. A MODEL HOST GALAXY: M33 

The nucleus of M33 contains a com pact star cluster that has 
elicited attention for a decade (e.g.. iKormendv & McClure 
U993llLauer et alJl998tfMatthews et al.ll999tlDavidgel2000l) 
because of its resemblance to the massive black hole inhab- 
ited nucleus of the Galaxy. No massive black hole, how- 
ever, has been found in M 33 (tMerritt. Ferrarese & Joseph! 
2001; Gebhardt et al. 2001). The interstellar medium and 
the magneti c field of M33 have also been extensively s tudied 
(e.g. IWilson & Scovillelll989t lEngargiola et"alll2003l iBeckl 
2000). Rotation curves of the more distant N SC host galaxies 
appea r to be rising even more slowly (|MaQh ews_& Gallagher 
120021) than that of M33 (ICorbelli & Saluccl2000l) 

We adopt the following para meters for M33: 
present NSC mass M cl ~ 1O 5 M (lLauer et alJ Il998t 
Long. Charles & Dubus 2002 ), surface d ensity of H I of 
E ~ IOM^pc"^ (lEngargiola et al J 120031). H I disk scale 
length /?o ~ lOkpc (Engargiola et al. 2003), average total 
magnetic field strength B ~ 6 /iG (Beck 2000), and a thermal 
sound speed of c, ~ 10 km s -1 . 

3. FORMATION OF NUCLEAR STAR CLUSTERS IN SITU 



^ Given a dynamical mechanism that transports gas from 

See http://www.ociw.edu/ociw/symposia/series/symposiuml/proceedings.html i ar g e r to smaller radii in the disk, an NSC can form from the 
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gas that accumulates at the center of the disk. The source of 
radial transport in non self-gravitating disks is the magnetic 
stress if the magnetorotational instability (MRI) develops in 
the disk. When there is differential rotation, torques arise be- 
tween concentric rings of interstellar gas. The Alfven velocity 
in a fluid with density p, velocity u relative to the circular orbit 
motion, and magnetic field B is given by = (47T/9) -1 /2 B. The 
few microgauss magnetic field measured in the disks of galax- 
ies such as M33 is compatible with being in equipartition 
with the observed velocity dispersion (~ 10kms~ ), where 
the density is estimated via p ~ Yjv/c s R (v is the velocity of the 
rotatio n curve of the galaxy at radius R). Sellwood & Balbus 
( 1999) pointed out that the velocity dispersion could be a re- 
sult of magnetohydrodynamic turbulence. 

Ignoring self-gravity of the fluid, the stress tensor Tr^ = 
(piiiRUfj, - uarUA(/>)) is an average of Reynolds and Maxwell 
stresses, where we have adopted cylindrical coordinates 
(R,(j),z)- Correlations between velocity and magnetic field 
components, which give rise to Reynolds and Maxwell 
stresses, induce an outward angular mom entum flux and an in- 
ward mass flux (Balbus & Hawlev 1998). The mass-weighted 
vertically averaged specific stress tensor is given by Wr^ = 
Tr<p/(p) and has dimensions of square velocity. 

We parameterize the average stress ten sor as (Wr^) = ac 2 , 
where a is a parameter introduced by Shakura & Sunvaev 
i ll 9731) that can be determined by measuring the stress-to- 
pressure ratio in numerical simulations. We assume that a 
and c s are wea k functions of radius. The ac cretion rate is then 
given by (e.g., Frank, King, & Rainel2002l) 
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47rac^ d 
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(R 2 X), 
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where k = [R~ 3 d / dR(R A Q 2 )] x l 2 is the epicycle frequency. 

The values of a associated with the MRI have 
been measured in three-dimensional simulations 
(e.g.. lHawlev. Gammie. & Balbusl Il995t IStone et all Il996t 
Armitage 1998). For the relevant isothermal equation of state, 
a ~ 0.1 when the initial (weak) magnetic field is uniform 
and perpendicular to the dis k plane and a ~ 0.01 when it 
is nonuniform. Recently, iKim. Ostriker. & Stond J2003I) 
conducted shearing-box simulations of a section of a gaseous 
galactic disk with an initial weak vertical magnetic field and 
measured a « 0.15-0.3. 

Atomic gas in a galactic disk can be treated as an ideal 
magnetic fluid if a sufficient fraction of ions are interspersed 
among the neutral atoms to provide for an inertial cou- 
pling of the neutrals to the magnetic field. iBlaes & B albus 
( 1994) showed that adequate coupling is achieved when 
the collision frequency of a neutral with the ions exceeds 
the local epicyclic frequency, 7/), > k, where 7 = (m n + 
m d~ l ( (7V )ni + C/pi i s expressed in terms of the density of ions 
Pi, the momentum rate coefficient (a v)„j « 1.9 x 10~ 9 cm 3 s -1 
( Draine, Roberge &Dalgarnoll"9 83). Here, m„ and m, are, re- 
spectively, the masses of the neutrals and the ions (taken here 
to be the hydrogen mass), and ( is the total ionization rate per 
molecule. 

IWolfireetaT1ll200l argue that in the Galactic disk, the ion- 
izing flux is made of extreme UV photons, soft X-rays, and 
cosmic rays, with a total of ( ~ 5 x 10~ 17 s" 1 . They also eval- 
uate the ionization fraction as a function of the density of neu- 
trals, x Jn) ~ 10~ 3 (»/10 cm" 3 )" 1 , consistent with earlier esti- 
mates (Kulkarni & Heilesl ll987l and references therein). In 
H I, the first term in the expression for 7 dominates. This 



implies that the MRI can develop if n < 2 x 10 s . Us- 
ing the available rotation curve data for M33 and borrowing 
the ionization fraction from the Galaxy, we estimate that k ~ 
10" 12 s" 1 atR=l pcandK-3x 10" 15 s" 1 atfl=l kpc. (Solid- 
body rotation would correspond to the maximum epicycle fre- 
quency of k ~ 3 x 10~ 12 s" 1 .) Therefore 7p, <C k throughout 
the disk, implying that ideal magnetic fluid is a good approx- 
imation to H I in disks of Scd galaxies, provided that condi- 
tions in these disks are similar to those in the Galactic disk. 

Our goal is to evaluate the gas mass accumulated at the 
center of the disk in time At. We assume that the total dy- 
namical mass enclosed within radius R in the center of the 
galaxy is a power-law of radius Md yn (R) = R{V 2 G~ l {R/R\y'~ 1 , 
where again 7 = -d\og p / d\ogr . We have parameterized the 
dynamical mass — including dark matter, stars, and gas — by 
the circular velocity Vi at some fixed radius Ri. Circular ve- 
locity at a generic radius is V(R) = V\{R / R\) x ~^ I 2 and thus 
fl, k oc RT 1 ! 2 . Surface density in the disk is in a quasi-steady 
state if the flux of mass through radius R is independent of 
radius, dM gRS /dR = 0. This requires that £ oc RT 1 ! 2 , and thus 
M gas = 2 7 rac 2 /?iE(/?i)yf 1 . 

The quasi-steady state is reached only at radii R < R dcc with 
accretion times f acc ~ 7r7? 2 S/M gas shorter than the age of the 
disk fdi s k, where R dcc is the maximum such radius. If we as- 
sume that the age of the disk of M33 is f^k ~ 1 Gyr and that 
the viscosity parameter is a ~ 0.1, we obtain R dcc ~ 500 pc. 

We compare the gas density profile S ce ntrai °c R~^l 2 just 
inside R RCC to the profile S exp oc e~ R l R ° obtained by ex- 
trapolating the exponential density profile that describes 
the disk at larger radii unaffected by accretion. Since 
(dT, centmi /dR)/(dT, exv /dR) = ~fR /2R. dcc at R acc , the density 
at R < R dcc exceeds the extrapolated density when 7 > 
2R dcc /Ro ~ 0.1. If the star formation rate (SFR) per unit 
gas mass is uniform in the disk, the excess gas density in the 
central few hundred parsecs implies an excess in the surface 
brightness profile of the stellar light. 

We proceed to evaluate the accreted mass. The gas mass 
accumulating at the center of the disk in time At equals 



M„ as - 10 7 M,r 



At 
10 9 yr 

( Rac 



0.1/ V 50 °P C 



lOkms" 

V(*acc) 



E(/?acc) 

10M Q pc" 



40kms 



(2) 



If ~ 10% of this gas is converted into stars, a cluster forms 
with a total initial stellar mass of ~ 10 6 M©. This is compati- 
ble with the mass of a typical NSC. 

Given that the development of MRI is contingent on differ- 
ential rotation in the disk, do disks in near solid-body rotation 
exhibit the same central accumulation? The instability can de- 
velop if the Alfven speed is not larger than the vertical sound 
speed times a function quantifying the degree of differential 
rotation (see Balbus & Hawlev 1998, eq. [110]). The condi- 
tion for instability can be written as ua < (2/Tr)q x l 2 c s , where 
q= \d\ogtl/d\ogR\. 

A disk could thus evade MRI if the rotation were solid 
body, i.e., if the total mass density were uniform. The growth 
time of the instability 2VT l q~ l , which here equals f grow ~ 
5 x 10 7 7- 1 (^/500pc)(V/40kms" 1 )- 1 yr, is short f grow < 
?disk even in very shallow (7 > 0.1) density profiles (note 
that p oc R' 1 ). Therefore the development of MRI and its 
consequences — the radially inward transport and the central 
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accumulation of gas — are universal in disks residing in sin- 
gular host dark matter halo profiles, however mild that singu- 
larity may be. After instability develops, the transport rate is 
independent of the dark halo density profile. 

4. NSC COULD NOT HAVE MIGRATED TO THE CENTER 

Here we evaluate the alternative scenario in which NSCs 
migrate to disk centers after forming elsewhere. We esti- 
mate the time for migration from radius R. Assuming that 
the host dark matter halo is spherical and concentric with 
the galactic disk, the total mass density inside the scale ra- 
dius of the halo is p(r) = poir/r^Y 1 + £oexp(-/?//?o)<$(z), 
where R = (jc 2 +y 2 ) 1 / 2 and 7, pp, r$, So, and Rq are con- 
stant parameters. Assume that the cluster of mass M c \ forms 
on a disk-planar circular orbit at distance R from the cen- 
ter. The timescale for orbital decay via dynamical friction 
is given by t ~ v^ irc /47rG 2 M / o<lnA ( Chandrase kharl 119431) . 
where v c ; rc is the circular velocity, p< denotes the density 
in dark matter particles moving slower than v c j rc , and In A 
is the Coulomb logarithm. The ratio p< j p is unknown as it 
depends on the phase-space structure of the dark halo. To 
derive a lower bound on the decay time we optimistic ally set 
P<r = P = (3-7)v 2 irc /47rGr 2 . Furthermore, Velazauez & White 
( 1999) found that In A ~ 2 for satellites sinking in galaxy disks 
with 7=1 halos. The decay time scale is then given by 



t > 3 x 10 9 yr 
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1 kpc 



( M d 
V 1O 6 M 



(3) 



Dynamical friction due to the disk is a more delicate is- 
sue. Beca use of the linear response o f disks to perturba- 
tions (e.g., iGoldreich & Tremaindll979T) . positive and neg- 
ative torques simultaneously act on the satellite. Negative 
torque is e xerted by the horseshoe orbits in the disk and 
amounts to (Ouinn & Goodman 1986) 

T^2\AB\H 4 R dZ/dR, (4) 

where A = -dfl/dlogR and B = A - ft are the Oort param- 
eters, and H = (GM c i/4A 2 ) 173 is the radius of the Hill sphere 
centered on the cluster within which the cluster potential dom- 
inates differential flow in the disk. 

Equation (0} is valid when the Hill radius is larger than the 
disk scale height h. For a smaller M c \ with H < h, the torque 
is smaller by fraction (H/h) 2 . If the rotation curve is "flat" 
fi oc R~ l , this result coincides with the detailed contribution 
from c orotation resonances evaluated by Donner & Sundelius 
( 119931) . Higher order resonances also contribute negative 
torque terms w hen R < Rp but their co ntribution is compa- 
rable or smaller ( Wahde, Donner & Sundelius 1996) than that 
of the corotation resonances and can be ignored. 

The inspiral time due to the torque in equation (0} is 

5/3 / \5/3 / R x-2/3 

^ , ,-, I!) ' 1 / > : I > CUV \ I iv 
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(5) 



where 7 henceforth denotes the negative logarithmic slope 
of the combined dark matter and baryonic density profiles. 
Deceptively, the timescale in equation (|5} formally becomes 
infinitesimal when the disk is in near solid rotation 7 — > 0. 
The implicit assumption H <C R breaks down in this limit, 
and we may try to replace H with R in equation @ to ob- 
tain / ~ td yR M c \/ 'MgasCR), where fd yn = /?/ v circ is the dynamical 



time and M gas is the enclosed gas mass. However, at radii 
where M c \/M g - ds (R) ~ 1, the outward horseshoe orbits that 
provide the negative torque are depleted, and the net torque 
thus changes sign, thereby reversing the sense of the migra- 
tion. 

We contrast the exceptionally long timescales derived 
above with the estimated stellar ages of NSCs. Using the 
bolometric luminosities and the density of stars on the asymp- 
totic g iant branch in the nucleus of M33, Stephens & Frogel 
(2002) inferred twin starbursts at ~ 0.5 and 2 Gyr. These 
ages could be reconciled with dynamical friction facilitated 
infall only by starting with very small initial ra dii (<C 1 kpc). 
Furthe rmo re, the best-fitting starburst mo dels of Gordo net al.l 
( 1999) and Long, Charles & Dubus (2002) imply the presence 
of another, very recent (~ 50-70 Myr) star formation episode 
yielding a stellar mas s of ~ 10 4 M,n or about 12% of the 
total mass. Similarly, B oker. van der Marel & VaccU ( 119991) 
determined the mean age of 10 6 8 — 10 78 yr for the nucleus 
of the giant Scd spiral IC 3 42 hosting a 6 x 1O 6 M cluster. 
Davidge & Courteau (2002) infer enhanced star formation in 
the recent 0.1 Gyr in the nuclei of nearby Sc galaxies NGC 
247 and NGC 2403. This evidence suggests that some NSCs 
are much younger (or contain a much younger component) 
than the migration times, and we thus reject the migratory hy- 
pothesis. 

5. DISCUSSION AND CONCLUSIONS 

When present, nonaxisymmetric features (bars and spi- 
ral structure) channel gas toward th e center of the galaxy 
dShlosman. Frank &~B egelmanl ll989l) . However, some NSC 
host galaxies exhibit negligible departures from axisymmetry. 
Self-gravity in a disk gives rise to an effective kinematic vis- 
cosity (e.g., iLin & Pringldll98"7l) . Stability of disk s against 
self-g ravity requires Q > 1 where Q = c s k/-kGY1 (fToomre 
1964). The disks of NSC host galaxies are largely sta- 
ble. Although that o f M33 contains approxim ately equal 
masses in stars and gas, Corbelli & Salucci ( 2000) infer Q > 1 
throughout, while our model implies Q 3> 1. Failure of the 
Toomre criterion to expla i n the SFR in M33 is well known 
(Wilso n & Scovilldll989i iMartin & Kennicuttl 12001 . Fur- 
thermore, Engargiolaet al. (2003) argue that giant molecular 
clouds (GMCs) in M33 are transients forming directly from 
the atomic gas with a characteristic mass ~7x 1O 4 M , which 
is much smaller than the Jeans mass c 4 /G 2 S ~ 3 x 1O 7 M 
of the disk. Therefore the GMCs are not like projectiles 
that facilitate radial transport through cloud-cloud interactions 
(e.g. JSilk & NormarJl981l) . 

We find that the radial transport of interstellar gas in galac- 
tic disks susceptible to MRI is sufficient to facilitate an ac- 
cumulation of gas in the central few hundred parsecs and the 
formation of a compact stellar system at the disk center. The 
sizes and the ages of the central stellar systems in this model 
correspond closely to those of the NSCs in Scd galaxies. 

We speculate that after a cluster forms in the disk, the sup- 
ply of gas flowing toward the nascent cluster is temporarily in- 
terrupted by the stellar winds and the supernova blast waves. 
During the life time of O and B stars ~ 5 x 10 7 yr, the stellar 
energy output fuels a hot low-density "superbubble" centered 
on the cluster. The radius o f the bubble reaches about t he disk 
scale height ~ 100 pc Ce.g.. lMac Low & McCravll988l) . After 
the stars with masses > 8M have been expended, fresh gas 
returns to the bubble and accumulates toward a subsequent 
starburst. The star formation history of an NSC is therefore 
punctuated by periodic starbursts separated by > 10 8 yr. 
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The radial transport of neutral gas in the disk implies a sur- 
face density enhancement inside a radius R icc measuring hun- 
dreds of parsecs. Remarkably, a "central light excess" a few 
hundred parsecs in radius, attributable neither to the NSC nor 
to the exponential stellar disk, has been observe d in severa l 
host galaxies (Boker et al. 2003) including M33 llKentll987l) . 
The excess is a natural consequence of the transport. Indeed, 
while the excess is observed in the stellar light and it is the 
gas that is transported by the MRI, the SFR in the overdense 
region may well trace gas density as elsewhere in the disk (Q 
is independent of radius for R < R acc ; see §[3}. 

Bright nuclear disk components of the late-type spirals de- 
void of luminous NSCs, such as in the superthin galaxy UGC 
7321 ({Matthews 2000), could also be the result of radial trans- 
port of interstellar gas. We expect that the NSC luminosi- 



ties depend strongly on the time elapsed since the most re- 
cent starburst episode. As dynamical data on NSCs become 
available, it will be possible to estimate the period between, 
and the magnitude of the consecutive episodes, which may 
lend insight into t he assembly of bulges and "pseudobulges" 
JKormendvll993l) . 
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